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THE RANDOM CONTINUED FRACTION TRANSFORMATION 


CHARLENE KALLE, TOM KEMPTON, AND EVGENY VERBITSKIY 


Abstract. We introduce a random dynamical system related to continued fraction ex¬ 
pansions. It uses random combination of the Gauss map and the Renyi (or backwards) 
continued fraction map. We explore the continued fraction expansions that this system 
produces as well as the dynamical properties of the system. 


1. Introduction 

In 1913 ( [Per50] ) Perron described an algorithm producing finite and infinite continued 
fraction expansions of real numbers of the form 

X = do-\ -, 

, ei 

di H- 

, . Cn-l 

d2+ • • +- 

dn + ■ • 

where do E Is and for each n > 1, e^-i G { — 1,1}, dn E N and dn + > 1. Moreover, in 

case the continued fraction is infinite, the algorithm guarantees that dn + en+i > 1 infinitely 
often. Perron called these expansions semi-regular continued fractions. 

Within this framework one can see a number of more familiar systems of continued frac¬ 
tions, each of which can be studied by a corresponding dynamical system. Regular con¬ 
tinued fractions, which correspond to letting each e„, = 1, are generated by the Gauss 
map Tx = ^ (mod 1). Backwards continued fractions, which were introduced by Renyi 
in |Ren57] . correspond to e„ = —1. Odd and even continued fractions ( |HKn2] ) and a- 
continued fractions ( |Nak81j ) can also be seen within this framework. In this article we 
define a random dynamical system which allows one to generate all semi-regular contin¬ 
ued fraction expansions of the form ([1]) for any given x and to study their dynamical and 
ergodic properties. We do not require the condition that dn + > 1 infinitely often, 

which makes our set-up slightly more general. 
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Over the last decade there has been a great deal of work on random dynamical systems. 
A random system is given by a hnite family of maps dehned on the same state space and a 
probabilistic regime for choosing one of these maps at each time step. The stndy of condi¬ 
tions that gnarantee the existence of an invariant measnre for snch systems was initiated 
by Morita f |Mor85] ) and Pelikan { |Pel84 ] ). They consider the case where each map in the 
random system is piecewise smooth with respect to some hnite partition and expanding 
on average. In |GB03] and |BG05] these resnlts are extended to when the probabilistic 
regime is position dependent. These results were further generalised by Inoue f |Inol2j i 
to more general underlying partitions, including ones with countably many elements. In 
[ANVj limit theorems are studied for random systems consisting of countably many maps. 
Their examples include families of maps that are piecewise smooth with respect to some 
hnite partition and are expanding on average. Rousseau and Todd studied hitting time 
statistics for random maps |RT15j . 


Random dynamical systems have also been used in relation to representations of numbers, 
see for example [Mor j . In particular, the introduction of the random /^-transformation 
has opened up new approaches to studying the dynamical and ergodic properties of /3- 
expansions and Bernoulli convolutions, see |DdV051 IDK131 IKeml3j . A key part of this 
approach was to study the invariant measures and ergodic properties of the random jS- 
transformation, as was done in |DdV05| IDdV07| IKeml4] . 


In this article we hrst introduce the random continued fraction map K. We show that 
it generates convergent continued fraction expansions for all points in its domain and we 
explore some of the properties of these expansions. The map K, and a related map R 
which is easier to analyse, present several challenges which are interesting from a purely 
dynamical point of view. In particular, R has countably many discontinuities and one of 
the interval maps dehning R has an indifferent hxed point. To overcome these difhcultieswe 
build on the work of Inoue |Inol2] . who studied transfer operators for countably branched 
skew product systems that are expanding on average. We can then use the results from 
|ANVj to obtain limit theorems. Other specihc examples of random intermittent systems 
have been recently considered in |BBD14j . Here Bahsoun, Bose and Duan studied limit 
theorems and mixing rates for random combinations of maps that are variations of the 
Manneville-Pomeau map. 


The paper is outlined as follows. In Section 2 we hrst show that every expansion of the 
form ([1]) is generated by K, and study further dynamical properties of random continued 
fractions. In Section 3 we prove that there exists an absolutely continuous invariant mea¬ 
sure for R. The density of this measure is of bounded variation. A key part of our approach 
here is to study the transfer operator associated with R. We show that R satishes condi¬ 
tions of Inoue |Inol2] which gives that the transfer operator is quasi-compact. In Section 
4 we use this to show that the invariant measure is fully supported. We can then employ 
results from |ANV] to obtain that R is mixing and that the Gentral Limit Theorem and 
Large Deviation Principle hold. Numerical evidence seems to suggest that the density in 
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fact is quite smooth. In the last section we discuss this in more detail and also mention 
some open questions and future directions. 


2. The random map 


2.1. Definition of the map. It is clear that x G M has an expansion of the form 


( 1 ) 


X = 


eo 


di + 




(^2 + ■ • + 


^n—1 


dn + 


where e^-i G {—1,1}, dn EN and dn + en > I for all n G N if and only if x G [—1,1]\{0}. 
If |x| > 1, we can subtract a suitable integer do from x and use the representation from 
m for X — do to obtain a continued fraction expansion of x. 


To hnd the right dehnition of the random dynamical system K, we hrst ask which x G 
[—1,1] have expansions that begin with a given choice of Cq and di. Let Cq and di > 1 
be given. Then x can be written in the form ([T]) if and only if x = for some y G 
[—1,1] \ {0}. This can be satisfied if and only if cq = sgn{x) and 

This typically gives two choices of the digit di. We see that 

^0 , 1 , 
y= -di = — - di. 

X |x| 

Similarly, if di = 1 then we require ei = 1 and have x = for some y G (0,1]. This can 
be satisfied for |x| G (|, 1] and we have 



As is standard with dynamical constructions of expansions of real numbers, the possible 
values of ei, d 2 are equal to the values of eo, di associated with y = -j^ — di. Thus we can 
generate all expansions of the form ([T]) using the following transformation. 

Let £ {0,1}} = {0,1}^ and let cr : —)■ be the left shift. We 

dehne the random continued fraction map K : x [—1,1] —)■ x [—1,1] by setting 
K{uj, 0) = (cr(a;), 0) and for |x| G (^, i], 
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Figure 1. The map Tq in red and Ti in blue. 


One can think of the map K as follows. Let maps Tq,Ti : [0,1] —)■ [0,1] be the Gauss and 
Renyi maps respectively, given by 


0 , 


if a: = 0, 


0 , 


( 2 ) Tqx = 


and Tix = 


— (mod 1), otherwise, 


X 


1 — X 


(mod 1), 


Let TT : O X [—1,1] [—1,1] be given by 7r{oj, x) = x. Then 

Tqx — cji, if X > 0, 


7r(K(LJ, x)) = 


Ti(x + 1) — coi, if X < 0. 


if X = 1, 
otherwise. 


Set Wo = 0 if X > 0 and 1 otherwise. This gives 

(3) 7r(Rr(w,x)) = T^o(x + Wo) - wi. 

Note that 7r(iL(w,x)) G [0,1] if Wi = 0 and 7r(iL(w,x)) G [—1,0] if Wi = 1. By iterating 
we see that for each n > 1 and all (w,x) G G x [—1,1], such that 7r(iF™’ (w,x)) 7 ^ 0 for 
0 < m < n, we have 


7r(K’^(u, x)) = (T^„_i o • • • o ZvJ(x + Wo) - w„ 


2.2. Random continued fraction expansions. For n = 1, set 


di = di(u, x) 


k + cji, if |x| G 


(— - 


CO 


if X = 0 
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and for n > 2, define dn{u;,x) = di(^K"’ We have 

(4) 7r{K{u,x)) = - di. 

and for each n > 1 snch that 'n'{K^{uj, x)) ^ 0 for all 0 < m < n, 


di + 7r(i^(a;, x)) 


, , (-!)"■ ~ 

^ 6^2 + 7r(i^(a;, x)) 


(- 1 ) 




di + 


(- 1 ) 


Ui 


(^2 + ■ • + 


(- 1 )^ 


dn + n{K^{uj,x)) 


The digit seqnence [dn{u,x))^^^ represents the continned fraction representation of the 
pair (cj, x) as given by K. If there is a smallest integer n snch that dn{oj,x) = oo, then 

{u, x) has finite random continued fraction expansion 


x = 


di + 


—(ZTjw 


(^2 + ■ • + 


d-n—l 


Now suppose that dn{u, x) is hnite for all n > 1. We want to show that (w, x) has infinite 
random continued fraction expansion 


X = 


- 1 ) 


^0 


di + 


- 1 ) 




do ■ ■ + 


- 1 ) 


^n — 1 


dn + ■ • 


For each n > 1, let — denote the convergents of the random continned fraction of (oj,x), 
i.e., write 


Qn 


di + 


(- 1 )“' 


(^2 + ' ■ ■ + 


dn 


2.3. Convergence of Partial Sums. To show convergence, we follow the approach for 
regular continued fractions with the important difference that in our case the numbers qn 
do not necessarily increase. As for regular continued fractions (see for example [DKn2] ). 
one can obtain the following relations: 

P-l = 1, Po = 0, Pn = dnPn-l + 

q-l =0, qo = 1, qn = dn^n-l + {-^T"~"qn- 2 - 
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Using these recurrences, induction easily gives that 

_ Pn+ 

Qn + qn-l'JT{K^iu),x)) 

Pn-lQn - Pnqn-1 = (-l)-(-l)- 0 +-+-n-l _ 

The next lemmas are needed to show that although the sequence {qn}n>i is not necessarily 
increasing, we still have lim^^oo ^ = 0- 

Lemma 2.1. For each n > 2, > 0 and if qn < qn-i, then Un-i = 1 and dn = 1, so 


Proof. We prove this by induction. For n = 1 we have qi = di > 1 = go- For n = 2, we 
have g 2 = d 2 di + (— 1 )'^^ and hence 

q2<qi ^ da < 1--—. 

di 

This can happen only if wi = 1, thus di > 1, and da = 1. This gives ga = di — 1 > 0 and 
the lemma. Now suppose the statements hold for all k < n, i.e., gfc > 0 and if g^ < qk-i, 
then Uk-i = 1 and d^ = 1, so = 0. If qn > 9n-l, then automatically qn > 0. So, suppose 
qn < qn-i- We have 

g„ = d„g„_i + (-ir"-ig„_a < g„-i ^ < 1 - 

Qn —1 

Since g„_a,gn-i > 0, g^ < g„_i implies that Un-i = 1 and hence d„_i > 1. The induction 
hypothesis then gives that < 1 and hence 1 < d„ < 2. This gives the lemma. □ 

Lemma 2.2. If qn < qn-i, then qn-2 < qn-i < qn+i- This implies that qn 7 ^ gn-i- 

Proof. By the previous lemma we have Un-i = 1, so g„_i > qn-2- We also have d^ = 1 
and thus Un = 0. Since qn > 0, this implies that g„+i = dn+ig^ + g^-i > gn-i- For the 
second part, note that if g„ = gn-i, then d„ = 1 — £ 2 ^ and by the previous lemma 

g„_i > qn- 2 - Since dn is an integer, this is not possible. □ 

Before we can prove that the process converges, we need a lower bound on the g„’s in case 
qn < qn-i- This is done in the next lemma. 

Lemma 2.3. Suppose qn < qn-i, so dn-i > 1. 

(i) Ifdn-i > 2 , then qn > qn-2- 

(ii) If dn-i = 2 and {uk, dk) = (1, 2 ) for all 1 < k < n — 1, then gn = 1 and g„_i = n. 
(hi) Suppose d„_i = 2 and there is al<k<n — 1 such that {uk, dk) 7 ^ (1, 2 ). Let k he 
the largest such index. Then qn > qk-i- 


(5) 

and that 

( 6 ) 
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Proof. Recall that Qn < Qn-i implies that hji = 1, = 0 and oon-i = 1, so dn-i > 1. 

Hence, = q^-i - qn- 2 - 

(i) If dn-i > 2, then 

Qn = dn-iqn-2 + ~ Qn-2 = (hn-1 — l)?n-2 + 

> 2 g „_2 + ( — 

If qn -2 > Qn-S, then qn > 2 g „_2 - qn -3 > qn- 2 - If qn -2 < Qn-S, then Lemma o gives that 
Un -2 = 0 and hence qn > 2 qn -2 + Qn-s > Qn-s > qn- 2 - 

For both (ii) and (iii), note that if = 2 and Uk = 1 for some 1 < k < n — 1, then 
Qk = 2qk-i - qk- 2 , so 

(7) Qk - Qk-i = Qk-i - qk- 2 - 

(ii) From (j7]) it follows that 

Qn Qn-l qn—2 d 2 Ql 2 ■ 2 -|- ( 1) 2 1. 

Moreover, for each l</c<? 7 , — lit holds that 

Qk = 2qk-i - qk-2 = Qk-i + (dfc-i - qk- 2 ) = Qk-i + 1- 
Hence, qn-i = n — 2 + qi = n. 

(iii) Let k be as given in the lemma, so {uk,dk) 7 ^ (1,2) and {uj,dj) = (1,2) for A: + 1 < 
j < n — 1. Then, 

Qn = Qn-i - qn -2 = Qk+i - Qk = dk+iQk + {-l)‘^'^qk-i -qk = qk + {-l)‘^^qk-i- 

If Uk = 0, then Qn = qk + Qk-i > Qk-i- If (-^k = 1, then dk > 3 and qk > Qk-i by Lemma [2Tl 
This gives 

qn = qk- Qk-i = {dk - l)qk-i + (-l)‘^'=-'gfc -2 > 2 gfc-i + (-l)‘^''"^dfc- 2 - 

As in the proof of part (i) we now get that if qk-i > qk- 2 , then qn > qk-i and if qk-i < qk- 2 , 
then Uk-i = 0 and we also get qn > Qk-i- D 

Proposition 2.1. Let x G [—1,1]\Q. For each oj, the digits dn{oJ,x) give a continued 
fraction expansion of x. 

Proof. For all n > 1 we have nsing ([5]) and ([ 6 ]) 



qn{pn +Pn-i'n-{K^{oj,x))) - Pn{qn + gn-ivr(A'”(a;, x))) 

Qn 

qn{qn + qn-lTT{K^{l-0,x))) 


7T{K'^{u,x)){qnPn-l - PnQn-l) 



qn{qn + qn-lTT{K'^{uj,x))) 


< 

1 

1 1 


\qn{qn + qn-l'n:{K^{uJ,x)))\ qn\qn-qn-l\ Qn 


We now show that lim^^oo — = 0. Let £ > 0. By Lemma Y272\ there exists a snbseqnence 

Qn 

{qnk)k>o snch that lim^^oo — = 0 and hence there exists an Ni snch that — < £. If there 
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is no n > A^i, such that Qn < Qn-i-, then we are done. If there is, let k be the smallest such 
index. If diVi+i) 7 ^ (1, 2), by Lemma 1^^ qi. > and then same holds for all other 

n > Ni. If dATi+i) = (1, 2), then set + 1. We then have ^ ^ < e. 

Moreover, for all n > we have > qk-i, since {uk,dk) = (0,1). This shows that the 
limit exists and is equal to 0. Hence we get 


for each x G [— 1 , 1 ]\Q. 


X = 


- 1 ) 


UJQ 


di + 


- 1 ) 




(^2 + 


- 1 ) 


UJ2 


da + ■ • 


□ 


3. Invariant densities 


Recall that Tq,Ti : [0,1] —)• [0,1] are the Gauss and Renyi map respectively. To study 
the dynamical properties of random continued fractions we replace the map K by the 
transformation R : fl x [0,1] —)■ G x [0,1] given by 

R{uj, x) = (acj, T^^x). 

The ergodic properties of R are easier to study and ([3]) gives the following simple relation 
between K and R. Recall that cuo = 0 if a; > 0 and 1 otherwise. Then, 

7 r ( R ”“^( a ;, Tqx)) = x)) + 

and hence, 

x)) =0JnR (-l)‘^"7r(i?”"^(a;, Tqx)) . 

Using this we can recover the digit sequence {dnioj, generated by K as follows. For 

(iv, x) G G X [0,1] define 

(8) b(io x)- I ^ e (aTT’ ’ 

\ cx), ifn;i + (-l)-ix = 0 . 

Write Ow for the sequence u' E fl satisfying oj'i = 0 and = ojn for all n > 1. Then for 

n > 1, 


(9) 


d„(a;, x) = 6 (R” ^( 0 a;,x)). 


Remark 3.1. The Gauss map Tq is intimately related to the well known Farey map, which 
is defined by 


X 


Fx = 


1 — X 
1 — X 


X 


, if 0 < X < 2 


2 ’ 


, if 2 < x < 1. 
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The Gauss map can be obtained from F by inducing on the hrst passage time to the 
interval [i, l]. For x G [0,1] and n > 1, set = 1 if x)) < | and = 0 if 

7r[K^~^{uj,x)) > If 'n'[R^~^{u,x)) = you can choose. Then n[R"'{u,x)) = F'^x for 
each n > 1. 


3.1. Existence. Let A denote the Lebesgue measure on [0,1]. We are going to show that 
R has invariant measures of type rup 0 ^p, where rup is the (p, 1 — p)-Bernoulli measure 
on G and Pp is a probability measure on [0,1] absolutely continuous with respect to A. 
For this we use |Inol2] . Recall that a function g : [0,1] [0,1] is said to have bounded 

variation if 

n 

y g ■■= sup ^ Ifir(Xi) - 5((Xi_i)| < CX). 

[0 1 ] 0=XO<Xl<X2<---<Xn=l 

The following is a very simplified, discrete version of the main theorem of |Inol2] : 


Theorem 3.1 (Inoue, |Inol2j L Given two non-singular maps Tq,Ti : [0,1] —)■ [0,1], let 
R : {0,1}^ X [0,1] —)■ {0,1}^ X [0,1] be given by 


R{u,x) = (cr(a;),T^i(x)). 


Let p G [0,1] and setpo = p andpi = 1—p. Fori G {0,1} let {R^k} be a countable partition 
of [0,1] into intervals and use int^R^k) to denote the interior of these intervals. Let g{i,x) 
be functions satisfying 


( 10 ) 


g{i,x) 


Pi 

TKx)\ 


on IJ^ int{R^k)- Assume that the following conditions are satisfied: 


(11) The restrictions of R to each interval int{R^k) are and monotone. 

(12) The weighted average expansion of R is uniformly positive for all x, i.e., 

sup ( 5 f( 0 ,x) + p(l,x)) < 1. 

3;e[0,l] 

(13) For each i G {0,1} the functions g{i,x) : [0,1] —)■ M are of bounded variation. 


Then there exists a probability measure Pp on [0,1] absolutely continuous with respect to 
the Lebesgue measure A with density function hp that is of bounded variation. Moreover, 
Up has the property that 

gp{A) = pPp{Tq^A) + (1 -p)/ip(TfM) 
for each Borel measurable set A C [0,1]. 


We have not included condition {A2) of Inoue since this automatically holds when one only 
has a finite number of maps R. We apply Theorem 13.II to our setting. 

Proposition 3.1. Suppose p G (0,1). Then the maps Tq and R satisfy the conditions 
(II), (12) and (13) of Theorem \3.1[ 
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Proof. For k>l, let /o,fc = and ■ Consider the following partitions 

of [0,1]: 

(11) 2io = {/o,fc}fc>i and Xi = {/i ,k} k>l- 

Then Tq and Ti are both and monotone on the interiors of the intervals of their respective 
partitions, so condition (II) is satisfied. Note that the fnnctions g{i, •) from (ITOH become 

g{0,x)=px‘^ and g{l, x) = (1 — p){l — xY. 

almost everywhere. We see that 

snp (^'(0, x) + 5 f(l, x)) = snp (x^ — 2(1 — p)x + (1 — p)) = max{l — p,p}. 

a:e[0,l] 3;e[0,l] 

So condition (12) is satished for all p G (0,1). Since both p(0, •) and g{l, •) are monotone 
fnnctions on the interval [0,1], we have 

V 9(0, ■)=P and \/ g{l, ■) = 1 - p. 

[0,1] [0,1] 

This gives (13), and the proof □ 

The conclnsions of Theorem 13.11 yield the following resnlt, the proof of which is standard, 
bnt inclnded for the convenience of the reader. 

Theorem 3.2. For any choice of parameter 0 < p < 1, there is an absolutely continuous 
probability measure Pp -C A such that the product measure mp ® Pp is invariant for R. The 
probability density function hp of the measure Pp is of bounded variation. 

Proof. Theorem 13.11 gives an absolntely continnons measnre Pp with the property that for 
each Borel set Ac. [0,1], we have 

( 12 ) pp{A) = pPp{Tq^A) + (1 - p)pp{T^^A). 

Take a cylinder [ji ■ ■ ■ jn] £ {0,1}^ and an interval (a, b) C [0,1]. Then 

■ ■ ■ Jn] X (a, b)) = [Oji ■ ■ ■ jn] X ((a, b)) U [Iji ■ ■ ■ jn] X Tf ^ ((a, b)). 

Hence, 

{mp^ Pp)(^R~^{[ji---jn] X (a, 6))j 

= Pmp{[ji ■ ■ ■jn])Pp(To\i(^, f»))) + (1 - P) "ip(bi • • ■jn])Pp(Ti\ia, 6))) 

= mp{[ji---jn])Pp{ia,b)) = (mp® pp){[ji---jn] x (a, &)). 

This gives the hrst part of the resnlt. The density hp is given by Theorem 13.11 □ 

If p = 1, then R rednces to the Ganss map and /ii(x) = If, on the other hand, we 

take p = 0, then R rednces to the map Ti, which has no absolntely continnons invariant 
probability measnre, bnt does have an infinite and cr-finite absolntely continnons invariant 
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measure with density /io(x) = K This was proved by Renyi in [Ren57j . The fact that R is 
expanding on average causes that hp is a probability density for all 0 < p < 1. It would be 
interesting to analyse the behaviour of hp as p —?• 0. 


3.2. Properties of the invariant measure. In this section we list some of the properties 
of the invariant measure rup ® ^p. In jInol2] Inoue proved the existence of the density 
hp by analysing a random version of the Perron-Frobenius operator. In our case, the 
corresponding operator £p : L^(A) —T^(A) is given by 

(13) {Cpf){x) = ^ ^ (h + ■ 

k>l 

Theorem 13.11 is proved in [Inol2] by showing that this operator has a hxed point in the 
space of functions of bounded variation, which is our function hp. Recall that functions 
of bounded variation can be modihed on a countable number of points to obtain a lower 
semi-continuous function. From now on we assume that hp is lower semi-continuous. Under 
the conditions of Theorem 13.11 the operator Cp is quasi-compact and constrictive. In this 
section we use these results to derive some dynamical properties of R. We hrst prove that 
R satishes a strong Random Covering Property. 

Proposition 3.2. Let I C [0,1] he a non-trivial interval. Then for every u ^ Ll, there is 
an n > 1 such that 

(T^„o---or^J/= [0,1). 


Proof. Recall the dehnition of the partitions Xq and Xi in the proof of Proposition 13.11 
If a non-trivial interval J is contained in one of the intervals in Xj, then RJ is again an 
interval and \{TiJ) > A(J). Hence, there is an m > 1, such that (X^^ o • • • contains 

an endpoint of one of the intervals in Tuim+i- This means that o o • • • o T^f)I 

contains an interval of the form [0, c) and an interval of the form (1 — c, 1) for some c > 0. 
Thus (X ^^^2 0 ---oT^J/= [0,1). □ 


It now follows that the measure /ip is in fact equivalent to A. The proof below is essentially 
the one from |ANVj . 

Proposition 3.3. Let hp he the prohahility density function from Theorem \3.S[ Then 
hp{x) > 0 for all x G [0,1). 


Proof. We know that hp is a function of bounded variation satisfying hp > 0, hpdpp = 1 
and Cphp = hp. Therefore there exists a non-trivial interval / C [0,1] and an a > 0, such 
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that hp > all- Fix a sequence a; = (a;i,a; 2 ,...) G ^2. Then for any n > 1, 
hp{x) = Cphp{x) > aCpli{x) 

— PuJi' ■■ PuJn ^ / \ 


> 


\{T 

(a;i,...,tj„)eO" ye(r„„o---or„^)-l{x} 

PoJl * * ■ 


o • • • o 




E 


\(T- 


^Myy 


o---oF^i)'2/r 

By Proposition 13.21 there is an n > 1 such that for all x G [0,1), 
Hence, for all 0 < p < 1 and all x G [0,1), hp{x) > 0. 


□ 


This proposition leads to the following result, the proof of which uses a standard technique 
that can be found for example in [BG97j . 

Proposition 3.4. The probability measure hp is bounded from above and away from 0. 


Proof. The fact that hp is bounded from above follows since hp is of bounded variation. In 
the previous proposition, we have established that hp{x) > 0 for all x G [0,1). Since hp is 
lower semi-continuous, it takes its minimum on [0,1]. Therefore, it is enough to show that 
hp{l) > 0. Let £ > 0 be small and let k>l. We consider part of the inverse image of the 
interval (1 — £, 1) under Tq (we could just as well use Ti). Note that 


and that 


A 


k + Vk + l — e 
1 1 


—) C T„-‘(l - E, 1) 


Hence, 


rA 


k + 1 k+1 — e 


k + Vk + l — e/J (fc-I-1)(A:-|-1 — e) 


3^)) < -^((1 ~ 1)) < (^ + 


k + 1 k + 1 — e 


It then follows that 

XunhJx) = 
xtl 


> 


> 


lim 

e^O A(^(l 


hp{x)dx = lim 


Pp{{l-e,l)) 


lim 

S —^0 


lim 


^ A((l - £, 1)) 

PPp{To\X - 1 )) + (1 - p)Pp{Tr\^ - 1 )) 


(^ + im(7TT>7TWi)) 


P 


hr 


««(fc + TA((;^,;pb)) (fc+TVfc+l 


> 0 . 


The fact that hp is bounded from above follows since hp is a function of bounded variation 
on a closed and bounded interval. □ 
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The fact that the operator Cp is quasi-compact on the set of functions of bounded variation 
allows us to obtain a number of consequences from the lonescu-Tulcea and Marinescu 
Theorem, as is done in many similar situations. The reader is referred to [BalOn] for example 
for an outline of this approach for deterministic maps. The spectral decomposition (which 
is already given in [Inol2] ) together with Proposition 13.41 gives that 1 is a simple eigenvalue 
for that there are no other eigenvalues on the unit circle and that £”/ —)■ j fdXhp in 
L^(A). This means that the system R satisfies the conditions of |AN Vj . (Note that R is not 
contained in the class of random Lasota-Yorke systems discussed in Example 2.1 of |AN V] . 
since both Tq and Ti have inhnitely many branches.) From |ANV] we then immediately 
get exponential decay of correlations: 


Proposition 3.5 (see |ANV] Proposition 3.1). There exist constants C > 0 and p < 1, 
such that for all functions f of bounded variation and all g G 


'[ 0 , 1 ] 


Cpf ■ gdpp 


'[ 0 , 1 ] 


fdpp 


'[ 0 , 1 ] 


gdpr 


< 


BV 


In particular, R is mixing. A proof can be found in [Mor85] for example. 
Theorem 3.3. The random transformation R is mixing with respect to rup ® pp. 


Proof. Define the function 0 : D x [0,1] —>■ [0,1] by 4>{uj,x) = uji + (—l)^^!^. Define the 
cylinder sets of order 1 by 


[d] X A(a)j 


[d\ X 0 


X 


a -]- 1 a. 


In general the cylinders of order k are given by [di • • ■ dk] x A(ai,..., ak)di---dk, where 


A(ai,..., ak)di-dk = n (^<^.-1 ° ° ^ ) • 

j=i 1 1 

Write [d] x A(a)| for [di ■■■dk] x A(ai ,..., ak)di-dk^ These sets form a generating semi¬ 
algebra for the cr-algebra on D x [0,1]. A straightforward computation using the previous 
proposition now gives that for any two such cylinders [d] x A(a)j and [c] x A(6)g, 

lim (rUp O pp){R~'^{[c] x A(fe)Q n {[d] x A(a)j)) 


where the convergence is in L^. 


= (rUp (g) Pp){[c] X A{b)l) {mp ® Pp){[d] x A(a) j) , 

□ 


We also immediately obtain the Central Limit Theorem and the Large Deviation Principle. 
Let 0 be an observable of bounded variation with 4>dpp = 0. Define Xk : D x [0,1] —)■ R 
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by x) = o • • • o x) and Sn = Y^k=o The asymptotic variance cr^ is 

defined by 

cr^ = lim — / S‘^d{mp® ^p). 
n^oo n 

In |ANV] it is proved that this limit exists in Proposition 3.2. Let M. be the set of all 
Borel probability measnres on [0,1] that are absolntely continnons w.r.t. A with a density 
of bonnded variation. 

Theorem 3.4 (Central Limit Theorem, see |ANV] Theorem 3.5). Let u ^ Ai. Then the 
process converges in law to A/'(0, under the probability measure mp 0 u. 

Theorem 3.5 (Large Deviation Principle, see jANV] Theorem 3.6). Suppose cr^ > 0. Then 
there exists a non-negative rate function c, continuous, strictly convex, vanishing only at 
0, such that for every u E A4 and every sufficiently small e > 0, 

lim — log ({mp 0 n){Sn > ne)) = —c{e). 

n^oo ' 

In particnlar these resnlts hold for A and Pp. 


4. Further Properties of the continued fractions 


Here we make a few observations abont the continned fractions prodnced by K. 


4.1. Counting Expansions. It is clear that K prodnces all semi-regnlar continned frac¬ 
tions of points X G [—1,1]. Since typically each seqnence u E Ll characterises a nniqne 
continned fraction expansion, K will generate nnconntably many different continned frac¬ 
tion expansions for most x. 

Proposition 4.1. Each x E [—1,1]\Q has uncountably many different continued fraction 
expansions given by K. Each x E ([—1,1] nQ)\{0} has countably infinitely many. 


Proof. To prove this we first show that x G ([—1,1] fl Q)\{0} if and only if for each u E Lt 
there exists an n > 0 snch that tt{K'^(lo, x)) E {^, —| : k E M}. One direction is clear, 
since if x is irrational then 7v{K{u, x)) is irrational for any u E Ll. Now assnme that 
X E [—1,1] O Q\{0}. First note that 7i{K(u,x)) E { — 1,0,1} if and only if x G | : 
k > l}. Write x = ^ with ro,ri G Z and |ri| < |ro|. If |ro| = |ri|, then the proposition is 
obtained with n = 0. If not, then 


ti{K(oj, x)) 


(-1)'^° (-l)‘^°ro - ridi 



[- 1 , 1 ]- 


Here r 2 G Z and |r 2 | < |ri|. If |r 2 | = |ri|, then we have the proposition with n = 0. If not, 
then by continning in the same manner we obtain a seqnence of integers tq, ri,..., with 
|r„| < • • • < |ri| < |ro|. This shows that the process mnst terminate after at most vq steps. 
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Now suppose X G [—1,1]\Q. Then there is no a; G and no n > 0, such that x)) G 

I ± i > l}. Hence, each uj produces a unique continued fraction expansion for x. This 
gives the hrst part of the proposition. 

For the second part of the statement, let | be an arbitrary element of ([—1,1] H Q)\{0} 
with p G Z, g G N and 0 < |p| < q. Then for each a; G H there is an < q, such that 
7t(^K^{u, |)) G {^, —I}- If Wn = 1 for all n > N, then there is some k > 1 such that 

(^dn(^uj, ^ j j = (di,..., (In, A; + 1, 2, 2, 2,...). 

If there is a smallest integer M >2 such that uj^+m = 0, then 

= {di,...,dN,k + l, 2, ■^.,2,1). 

M—1 times 

Finally, if wat+i = 0, then 

(dn(ui,^^'^ = {di,...,dN,k). 

Hence, K generates only countably many different expansions for x. □ 


Proposition 14.11 gives all the possible endings of digit sequences of rational points | gen¬ 
erated by K. From Figure [T] it is clear however that | also has expansions that are not 
generated by K. The corresponding digit sequences are 


{^di,..., dj\f, k 1,1), (^di,..., dj\f, k 1, 2,2, 2,...), 


and for each M > 2, 


(di,..., dN, k — 



M—l times 


The fact that K does not generate these expansions is a consequence of having dehned K 
by taking the intervals (^;^, left-open and right-closed. So all points in [ — 1,1] fl Q 
have expansions of the form ([T]) that are not generated by K, but we know exactly what 
the missing expansions are. 


a-continued fractions. The a-continued fraction transformation was hrst introduced 
Nakada in |Nak81] . Given a G [0,1], dehne the map Tq, : [a — 1, a) ^ [a — 1, a) by T^O = 
and 


= 


X 


- (a - 1) 


for X 7 ^ 0. First note that for x < 0, Ti(l -f- x) = |i| — [|y|J. Also note that 


by 
= 0 


1 

_ 

1 

< a 

1 

— O! -f- 1 

— 

1 

X 


X 


X 



X - 


and 


1 


1 

> a ^ 

1 

— -|- 1 


1 

— 

— 

— 

— 


— 

X 


X 


X 



X - 
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We can define a sequence oj = (<Vn)n>i G such that for each n > 0, x)) = T^x 

as follows. Note that if 0 < x < a, then either Tqx G [a — 1, a) in which case we set oii = 0 
or Tqx — 1 G [a — 1, a), in which case we set oii = 1. Similarly for a — 1 < x < 0 and the 
map Ti(x + 1). Then 

r^o(x + Uq) -Ui = T„x. 

We can apply the same procedure for every n > 2 with the point o • • • o T^q){x + 

Wo) — ojn- One can easily show, using (Ej) and induction, that for each x G [a — 1, a) and 
n > 1, one gets tt[K^{u,x)) = T^x. 

4.2. Restrictions on digits. For regular continued fractions, expansions with restrictions 
on the digits have been thoroughly investigated. Continued fractions with bounded digits 
have applications in many helds, ranging from formal language theory to diophantine 
approximation to pseudo-random number generators, see |Sha92] for a survey on this topic. 
For random continued fractions the following is easily observed. 

Proposition 4.2. Let A G N be such that it does not miss two or more consecutive integers, 
i.e., if k ^ A, then k — 1, k + 1 E A. Then all numbers in [—1,1] have a continued fraction 
expansion of the form m with dn E A for each n. 

Proof. From Figured] it is clear that if we delete the branches \ — k and —^ — k for k E A, 
then all the other branches still cover all of the interval [—1,1]. This means that for each 
X G [—1,1] we can choose an w G such that dn{oJ, x) ^ k for all n > 1. □ 

From this proposition it follows immediately that each x G [—1,1] has a continued fraction 
expansion with only even or only odd digits. Moreover, note that if we remove all branches 
corresponding to the odd (or even) digits, then there is no overlap in the system except on 
the points ±^, implying that each irrational x G [—1,1] has a unique expansion with only 
odd (or even) digits. 

Regular continued fractions with bounded digits are very well studied. It has been known 
since the works of Jarnik f [Jar32] i and Good f [Goo41] i that for each N > 2 the set of 
points X G [0,1] that has a regular continued fraction expansion with digits not exceeding 
N has Hausdorff dimension strictly between 0 and 1. Many people have given estimates 
for the Hausdorff dimension of the set Ei ^2 of points with a regular continued fraction 
expansion using only digits 1 and 2, with the most recent result by Pollicott and Jenkinson 
f |JP01] ) calculated the hrst 25 digits of the Hausdorff dimension of this set: 

dimH(F^i,2) = 0.5312805062772051416244686.... 

In the random case we have the following. 

Proposition 4.3. The set of points x E [—1,1] that have a continued fraction expansion 
of the form (Qj) with dn E {1,2} for all n > 1, has positive Lebesgue measure. In fact, it is 
a countable union of intervals, containing the interval (^, l] ■ 
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Proof. Let A 12 denote the set of points that have an expansion of the form ([T]) with 
only digits 1 and 2. Then ^41^2 hi ( — |, |) =0. First note that the points ±^, k = 
2,3, have conntably many expansions nsing only 1 and 2 as discnssed after the proof of 
Proposition 14.11 These expansions are not generated by K. For all other points in Ali 2 
snch expansions are generated by K. We consider these expansions, see Figure |2] for an 
illustration. 


The hrst digit of any x G Ai ^2 hi ((|, U [ “ “!)) i^^st be 2, coming from uji = 0. For 

n > 1, let Sn denote the map SnX = |-| — n. Then for points in the set 


s=Uv(o. 


n>l 


n 


1 1 
3 ’ 2 


u 


1 

2 ’ 


K does not generate an expansion with only digits 1 and 2. Note that S 2 has hxed point 
\/2 — 1 and > 1 for all t G U [ — —|) • Since 


\/2 - 1 ^ 


“' 3 ' = 


3 1 

7’2, 


U 


1 _3 
2’~7 


and 




°’3)- 


2 7 


5’ 17 




7 2 


17’ 5 


it follows that E hi [ — and that is a countable union of intervals. For 

the set (|, 1 ] n [ — 1, — I) we have the following. Since 




-1 


o.b = u.i 


u 




and 




-,o 

3’ 


- 

rl 3\ / 3 1] 



. 

L 2 5/ V 5 2 J 


we see that on l] and [ — 1, —|) we will have to use S 2 and on (i, |) and ( — f, —^) 
we will have to use Si. Note that 


^ri7n(-,i 


c 


2 5 
L3’7, 


Since | < | < | < |, on all of the interval l] we can avoid being mapped into E by 
making appropriate choices between Si and 82 . Hence 


H-1,2 — 


L3 


,1 


U 


1, 


\E. 


□ 


In 1994 Lehner ( [Leh94] ) proved that every irrational number x in the interval [1,2) has a 
unique continued fraction expansion of the form 

X = 60 H-, 


‘ • 

where (e„,6„) equals either (1,1) or (—1,2). In [DKOO] Dajani and Kraaikamp established 
a relation between these Lehner expansions and the regular continued fraction expansions, 
produced by the Gauss map. Here we see that all points in Hi 2 have a Lehner-like expansion 
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-1 . 1.1 0 1 1 1 
2 3 3 2 


Figure 2. Points that have a continued fraction expansion generated by K 
with only digits 1 and 2 only use the maps x i—)■ |i| — 1 and x i—)■ |^| — 2. 


with bo = 0. Moreover, there are points that have more than 1 such expansions. This holds 
for example for all points in the interval (|, |) or in the interval |). 

In spite of the result from Proposition 14.31 we have the following results on the asymptotic 
geometric and arithmetic mean of the digits, which are similar for the regular continued 
fractions. 


Proposition 4.4. For {nip 0 iJ,p)-a.e. {oj,x) G x [0,1], we have 

1< lim (di(a;, x) • • • < oo, 

n—^oo ^ 

and 

di{u,x)-\ - \-dJoj,x) 

hm -= oo. 

n^oo 77 , 


Proof. We use the ergodicity of the map R. Recall the dehnition of the function b : 
12 X [0,1] —)■ N U {cxd} from ([H]). For the hrst statement, we hrst show that log 6 G L^. 
Suppose hp{x) < M for all x G [0,1]. Then, 


'Ox[0,l] 


\ogbd{mp Up) < / logndx + p{l — p) / log(n + 1) 


dx 


n>l 


+p{l - p) 


n+1 

n 

r n+1 
I n-1 


logndx + (1 — p)" 


n+1 

n 

r n+1 

I n-1 


log(n + 1) dx 


plogn ^ (1 — p) log(n + 1) 


n>l 


n{n + 1) 


n{n + 1) 


< oo. 
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There also is an m > 0, such that hp{x) > m for all x G [0,1]. Similarly as above, we 
obtain 


'Ox [0,1] 


log b d{mp ® Up) > m 


plogn ^ (1 — p) log(n + 1) 


n>l 


n{n + 1) 


n{n + 1) 


> 0 . 


Then by the Birkhoff Ergodic Theorem, for {rrip (8) pp)-a.e. {u, x), 


1 ? \ f 

0 < lim —log b(R\u,x)) = / 

Jn 


i=0 

The result now follows from 


Ox [0,1] 


logbd{mp (8) fip) < oo. 


For the other statement, one hrst notices that b{u, x) > i —1 if oii = 0 and b{u, x) > —1 

if oil = 1. Hence, b ^ L^. Therefore, consider the functions 

[ 1 i](x), ifn;i = 0, 

bN{uj,x) = < , / / \ -f 1 

I if wi = 1. 

Then each b^ is bounded and the sequence is increasing, so by Beppo-Levi’s Theorem, 


(14) 


lim 

N^OO 


bdijUp (8) Up) = oo. 


/ bn d{mp 0 pp) = 

/ox[o,i] ,yox[o,i] 

Moreover, the Birkhoff Ergodic Theorem gives the existence of a {nip 0 /ip)-measure 1 set 
of {uj,x), such that 

n-l „ 

bjM d{nip 0 /ip) 


1 ' ; r 

lim —'S^bN{R\uj,x)) = / 
1^00 n /o 

i=0 


Ox [0,1] 


for all > 1. Then by flTT)) we get that for {nip 0 pip)-SL.e. (ca, x). 


n—1 


lim inf — di{u,x) = liminf — 6(/2*(Oc<;, x)) 

n^oo n n^oo fi ^ ^ ^ 

2=1 2 = 0 


> 


1 ? \ f 

lim lim — N b]\f{R\0uj,x)) = / b^ d{nip® pp) = oo. 
N^oon^oon^^ iox[ 0 ,l] 


This gives the result. 


□ 


5. Further Questions and Comments 


So far we have established that the density hp is of bounded variation and is bounded away 
from zero. It satisfies the equation 


hp{x) /lphp{x) ^ ^ 

k=l 


P 


h I 

{k + x)'^ ^ \k + X J ' {k + x) 


+ 


P 




1 - 


X + k 


= pCohp{x) + (1 - p)Cihp{x), 
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where Cq, Ci are transfer operators for Tq and Ti respectively. Both operators Cq and Ci 
preserve cones of positive smooth (analytic) fnnctions on [0,1]. Considerable effort has 
been pnt in the analysis of the spectral properties of £o, which are now well nnderstood, 
see |Iosl4] for a recent overview. It would be interesting to study the spectral properties 
of Cp in more detail. Based on simulations, we suspect the following. 

Conjecture 5.1. For each 0 < p < 1 the function hp is a smooth function on [0,1]. 


The exact smoothness condition is to be determined. We presume that applying relatively 
standard techniques one could strengthen the results of the present paper by showing that 
hp G (^^([0,1]) for all /c G N. In fact, we believe that the density is (7°°; however, most 
probably it is not real-analytic (which is the case for the Gauss map). 


This conjecture is motivated as follows. In |MR87( IJGU03] the authors investigate the 
relation between Cq and the the integral operator /Cq acting on the Hilbert space L^(M+, p), 
given by 

^o0(s) = / dfi{t), 

Jo \ st 

where Ji is the Bessel function of the hrst kind, and p is the measure on M_|_ with the 
density 

7 ^ 7 

du = —- dt. 

e* — 1 

The operator /Cq has a symmetric kernel Kq{s, t) = , and has several nice properties, 

e.g., is nuclear. In a similar fashion, existence of a positive smooth hxed point of Cp will 
follow from the existence of a positive fixed point of /Cp 

ICp(j){s) = p/Co0(s) + (1 - p)/Ci0(s) 


(15) 


= P 


Ji(2\/sf) 

\/sf 


0(t) dp(f) + (1 - p) 


/i(2\/gf) 

\/sf 


(j){t)e *dfj,{t), 


where R is the modified Bessel function of the hrst kind. Technical difficulties arise from 
the fact that the kernel Ki = albeit monotonic and positive (c.f., Kq is oscillating), 

is not integrable. 


An exact formula for the density hp would of course settle the conjecture. The most 
successful approach to constructing invariant densities for continued fraction transfor¬ 
mations has been to build natural extensions of the transformations, see for example 
[NakSll IHaa021 IIS061IDKS091 IKSSlOl IKSS12j and the references therein. This approach 
was also effective in determining the invariant density of the random /5-transformation (see 
[Keml4] ). in which case the invariant density was not just a linear combination of the 
invariant densities of the two maps making up the random transformation. We have so far 
not been able to build a natural extension of the random continued fraction map. 


With an expression for the density one could study frequencies of digits for the continued 
fraction expansions. It would also be interesting to consider the size of subsets of [—1,1] 
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obtained like those in Propositions 14.21 and 14.31 For example, for any two consecutive digits 
n and n + 1 is there set with positive Lebesgue measure such that all points in this set 
have an expansion using only these digits? 

Even without a formula for the density, one could study the behaviour of hp as a function 
of p. Of particular interest is the case when p —)■ 0, since ho is unbounded and hp is of 
bounded variation for each p > 0. A similar question can be asked for the metric entropy 
of R. Can one calculate this? The Shannon-McMillan-Breiman Theorem with the cylinder 
sets from the proof of Theorem 13.31 can be of help here. How does it behave as p —)■ 0? 
This last question could be seen as an analogue of the question of studying the entropy of 
a continued fractions as was done in |CT13] . 
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